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ABSTRACT 

We present spatially-resolved dynamics for six strongly lensed star-forming galax- 
ies at z = 1.7-3.1, each enlarged by a linear magnification factor ~ x8. Using the 
Keck laser guide star AO system and the OSIRIS integral field unit spectrograph 
we resolve kinematic and morphological detail in our sample with an unprecedented 
fidelity, in some cases achieving spatial resolutions of ~100pc. With one exception 
our sources have diameters ranging from 1-7 kpc, integrated star formation rates of 
2 — 40 M Q yr _1 (uncorrected for extinction) and dynamical masses of io 9 - 7-10 - 3 M Q . 
With this exquisite resolution we find that four of the six galaxies display coherent 
velocity fields consistent with a simple rotating disk model. Our model fits imply ratios 
for the systemic to random motion, V c smi/<7, ranging from 0.5-1.3 and Toomre disk 
parameters Q < 1. The large fraction of well-ordered velocity fields in our sample is 
consistent with data analyzed for larger, more luminous sources at this redshift. We 
demonstrate that the apparent contradiction with earlier dynamical results published 
for unlensed compact sources arises from the considerably improved spatial resolution 
Os . and sampling uniquely provided by the combination of adaptive optics and strong 

gravitational lensing. Our high resolution data further reveal that all six galaxies con- 
tain multiple giant star-forming Hn regions whose resolved diameters are in the range 
300 pc - 1.0 kpc, consistent with the Jeans length expected in the case of dispersion 
/\ ' support. From the kinematic data we calculate that these regions have dynamical 

masses of l0 8 - 8 ~ 9 - 5 M Q , also in agreement with local data. However, the density of 
star formation in these regions is ^100x higher than observed in local spirals; such 
high values are only seen in the most luminous local starbursts. The global dynam- 
ics and demographics of star formation in these Hil regions suggest that vigorous 
star formation is primarily governed by gravitational instability in primitive rotating 
disks. The physical insight provided by the combination of adaptive optics and gravi- 
tational lensing suggests it will be highly valuable to locate many more strongly-lensed 
distant galaxies with high star formation rates before the era of the next-generation 
ground-based telescopes when such observations will become routine. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift 



1 INTRODUCTION 

Studies of star-forming galaxies at high redshift (z > 1.5) 
over the past decade have mapped the demographics of the 
populations as a whole yielding valuable data on their inte- 



grated star formation rates, stel lar and dynamical masses , 
metallicities and morp hologies ilShapley et alj l200ll. 120031 : 
iReddv fc Steidell 12004 : lErb et alj 120061 : lLaw et alj l2007al ). 
From these data, a consistent picture of the development 
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of the comoving density of star formation and its con- 
tribution to the pres e nt day stellar density h as emerged 
[jDickinson et al.ll2003l : iHopkins fc Beacomll2006h . Following 
a rapid rise in activity at early times, corresponding to the 
redshift range 2 < z < 6, the star formation rate has decline 
markedly over the past 8 Gyr corresponding to the interval 
< z < 1. Of particular interest is the relationship between 
sources observed at the time close to the peak epoch of star 
formation, 2 = 2 — 3 and the population of massive galaxies 
observed today. 

Recent theoretical work has focused on the mechanisms 
via which early galaxies assemble their stars and evolve 
morphologically to the present day Hubble sequence. Mi- 
nor mergers and cold stream gas accretion have emerged as 
a possible means of building disks, central bulges and ellip- 
tical galaxies from ear l y star forming Lym an break galax- 
ies i|Dekel et al.l 120091 ; ]Br"ooks et al.l l2009h . These models 
of galaxy evolution rely on simple descriptions of complex 
physical processes such as gas cooling, star formation and 
feedback mechanisms. It is therefore crucial to undertake 
relevant observations to constrain these processes. The most 
direct observational route to making progress is high quality 
resolved spectro-imaging of early galaxies which can be used 
to determine both their dynamical state and the distribution 
of star-formation. 

Integral field unit (IFU) spectroscopy of nebular emis- 
sion lines is now yielding valuable kinematic data for a grow- 
ing sample of luminous high-redshift star forming galax- 
ies llGenzel et aij|2006t 12008: JForster Schreiber et al.ll200r3 . 
120091 : lLaw et al.1 l2007bl . 120091 ) . Although the galaxies ob- 
served so far do not yet in any way constitute a complete 
well-defined sample, some important results have emerged. 
The data reveal a mix of dispersion-dominated systems, ro- 
tating systems, and major mergers with the common obser- 
vation that all galaxies studied so far show relatively high 
velocity dispersions of ~ 50 — 100 kms -1 . A larger frac- 
tion of the more massive (Mdyn > 10 11 Mq) galaxies ap- 
pear to be rotationally supported, with lower mass galax- 
ies (Mdyn ^ 10 10 Mp)) tendin g to have dispersion-dominated 
kinematics l|Law et al.ll2009l ). 

However, these studies are hampered by the poor spa- 
tial resolution inherent in studies of distant star-forming 
galaxies whose typical sizes are ~ 1-5 kpc. Even with adap- 
tive optics on a 8-10 meter aperture, the physical resolu- 
tion is limited to FWHM > 1 kpc which means only a 
few independent resolution elements. Consequently it is un- 
clear whether the velocity shear observed in the smaller 2- 
5 kpc sources arises from rotation or merging. Moreover, al- 
though there are suggestions that distant star-forming re- 
gions have sizes much larger than those of local Hn re- 
gions l|Elmegreen fc Elmegreer]|2005l ). the claim remains un- 
certain since the distant regions are not properly resolved. 
Thus, the diffraction limit of the current generation of op- 
tical/infrared telescopes limits progress, even with all the 
adaptive optics tools of the trade. 

Fortunately, sources which undergo strong gravitational 
lensing are enlarged as seen by the observer and can thus 
be stu died at much higher physical resolution in th eir source 
plane (jNesvadba et al. 2006; Swinbank et~ai1l2007h . This im- 



provement in physical resolution enables us to better dis- 
tinguish various forms of velocity field and to examine the 
properties of star-forming giant Hn regions. A convincing 
demonstration of the b enefits of le n sing a s applied to this 
topic was the work of IStark et al.l (|2008T ) which analyzed 
IFU data on a, z — 3.07 galaxy magnified in area by a factor 
of 28 x demonstrating resolved dynamics on ~100 pc scales. 
This represents an order of magnitude improved sampling 
compared to the earlier studies cited above. 

Here we present IFU observations of a further five 
strong ly lensed ga l axies which we add to the source stud- 
ied bv lStark et al.l {2008). In addition to studying the well- 
sampled kinematic structure of a population of sub-L* star- 
forming galaxies at z ~ 2, our fine resolution allows us to 
examine the size, dynamical mass, and luminosity of typi- 
cal star-forming regions, thereby probing the conditions in 
which most of the star formation takes place. We find that 
all of our galaxies contain multiple, distinct star forming re- 
gions. These results offer a preview of the science which will 
be possible with 30 meter class optical/near-IR telescopes 
with a diffraction limit comparable to the resolution of our 
data; they provide the most detailed view to date of typical 
high-redshift star forming galaxies. 

A plan of the paper follows. In §2 we discuss our se- 
lection of targets and Keck observations, including the pro- 
duction and use of gravitational lens modeling essential for 
achieving a high resolution in the source plane. In §3 we 
analyze our spectroscopic results, discussing each object in 
turn and summarizing the kinematic data. In §4 we examine 
the physical scale of our star forming regions in compari- 
son to those seen locally. We summarize our conclusions in 
§5. T hroughout we use a WMAP cosmology (|Spergel et al.l 
l2003h with ft A =0.73, fi m =0.27, and H =72kms" 1 Mpc" 1 . 
All quoted magnitudes are in the AB system unless other- 
wise noted. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target Selection 

Ideally, in considering the selection of high redshift galaxies 
for more detailed resolved studies, care would be taken to 
construct a mass-limited sample, perhaps defined addition- 
ally according to the integrated or specific star formation 
rate. In practice, even with 10 meter class telescopes, studies 
are limited by the surface brightness distributi on of gaseous 
emiss i on lines as discussed by earlier workers ([Cenzel et al 



I 2006L l2008l: JForster Schreiber et al.1 l200d . 120091 : lLaw etal 
l2007bl . l2009h . In examining gravitationally-lensed systems, 
which have the unique advantages of probing to less lumi- 
nous (and presumably more typical) systems, and offering 
an increased spatial resolution in the source plane, an addi- 
tional criterion is the reliability of the mass model used to 
invert the observed data into that of the source plane. 

Over the past decade, using the growing ar chive of HST 
ACS and WFPC-2 images of rich clusters (e.g. ISantos et al.l 
12004 ISand et al.ll2005l ; IStark et alj|2007l ) and through spe- 
cific projects such as the MAssive Cluster Survey (MACS, 
lEbeling et ai1l200lf ) and Local Cluster Substructure Survey 
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(LoCuSS, ISmith et alj|2005l ), we have studied a large sam- 
ple of lensing clusters and undertaken a search within these 
HST images for promising magnified galaxies (arcs) at high 
redshift. The clusters presented herein do not in any way 
represent a well-defined sample and, inevitably, there are 
biases to those which act as spectacular gravitational lenses. 

The distant sources described here were located and as- 
sessed via two sequential programs. Firstly, a comprehensive 
multi-object spectroscopic campaign involving Keck and the 
VLT was used to verify the redshift of the candidate arcs as 
well as to determine the redshifts of other multiply-imaged 
sources in e ach cluster in order to construct robust mass 
models (e.g. iRichard et al.l [20071 ') . As the follow-up spec- 
troscopy was undertaken at optical wavelengths, for those 
sources beyond z ~2, this provided valuable insight into the 
strength and distrib ution of Lya emissio n. However, as with 
previous work (e.g. IShaplev et alj|2003h we found that the 
Lya line did not necessarily give a good indication of the rel- 
ative strengths and extended nature of those nebular emis- 
sion lines ([On], H/3, [Oni], Ha) critical for resolved studies 
with laser guide-star assisted adaptive optics (LGSAO). Ac- 
cordingly, we used the Keck II spectrograph, NIRSPEC, in 
seeing limited mode in a second campaign to 'screen' all 
promising targets in advance of undertaking the more de- 
manding resolved studies with the LGSAO-fed integral field 
unit spectrograph OSIRIS. Once the redshift of the lensed 
source was known from optical spectroscopy, we usually fol- 
lowed up the source in more than one near-infrared line, in 
order to ascertain which would be the most economic tracer 
of the velocity field. 

About 75% of our initial sample of 30 z > 1.5 arcs have 
so far been screened and have suitably close reference stars 
for tip-tilt correction. In order to achieve an adequate signal 
to noise with OSIRIS, we used the NIRSPEC data to limit 
the sample to those with an emission line surface brightness 
> 10~ 16 ergs _1 cm _2 arcsec -2 in a region uncontaminated 
by night sky emission. As a guide, this surface brightness 
criterion corresponds to a SFR density > 0.3 M0yr _1 kpc -2 
for Ha at z = 2. Our sample of 6 sources is drawn from a 
list of ~ 20 objects fulfilling these criteria. A full account 
of the NIRSPEC screening program and a discussion of the 
dis tribution of line stre ngths found in that survey is given 
bv lRichard et all (|2009l ) . Key properties for the sample of 6 
object s, which includes the source discussed bv lStark et al.l 
(2008), are given in Table 1. The ACS images of our lensed 
targets are given in Figure [1] and we briefly discuss these 
below. 

C10024 + 1654 Deep optical and near- infrared imaging 
of the z = 0.39 cluster CI 0024+1654 with HST reveals 
five distinct images of a background galaxy originally iden- 
tified by its blue color and clumpy ring-like morphology 
l|Collev et al 1 ll996l . Figure [TJ . Optical spectroscopy estab- 
lished a redshift of z = 1.675 from interstellar absorption 
lines ( Broadhu rst et all |2000| ) . Near-infrared spectroscopy 
shows bright, spatially extended Ha emission at z = 1.6795 
(Richard et al. 2009) . To map the velocity field, we observed 
the counter image northwest of the cluster with OSIRIS 
around the redshifted Ha and [Nn]A6583 emission. 

MACS0451+0006 The HST ACS V- and /-band imag- 



ing of MACSJ0451+0006 shows a striking, multiply-imaged 
arc ~30 arcseconds east of the brightest cluster galaxy 
(Figure [1]). Follow-up spectroscopy with the FORS multi- 
object spectrograph on the VLT confirmed a redshift of 
2 = 2.008 via the identification of Lya and UV ISM absorp- 
tion lines (SinA1260, OiA1303, CivA1549 and SiivA1400) 
and HellA1640 and Cm]A1908.7 in emission (PID: 078.A- 
0420). Follow-up observations targeting the nebular Ha 
emission with the NIRSPEC longslit spectrograph on Keck 
II showed stro ng, spatially extend ed Ha at a redshift of 
2.0139±0.0001 (|Richard et alj|2009h . 

MACSJ0712+5932 HST V- and /-band imaging of 
MACS J0712+5932 reveals a prominent triply-imaged back- 
ground galaxy. Each of the galaxy images clearly com- 
prises two prominent clumps with at least one fainter re- 
gion with surface brightness fj,y > 23 mags arcsec -2 within 
a diffuse halo (Figu res 1 1 1 2 1) . Follow-up near-infrared spec- 
troscopy l|Richard et al.l l2009) shows strong nebular Ha, H/3, 
and [Oni]AA5007, 4959 emission at z = 2.6462, concentrated 
in the UV-bright regions. 

MACS J0744+3927 HST V- and /-band imaging of 
MACS J0744+3927 reveals an arc approximately 15" west 
of the cluster center. The optical redshift z = 2.207 ob- 
tained via identification of strong Lya absorption and mul- 
tiple rest-frame UV absorption lines was later confirmed by 
the detection of strong Ha and [Nil] emission at z = 2.209 
with NIRSPEC. 

C10949+5153 The HST U-band image of CI 0949+5153 
shows a clumpy, elongated arc resolved into two diffuse com- 
ponents (Figure [TJ. Further spectroscopy around the red- 
shifted nebular emission lines shows extended Ha, H/3 and 
[Om]AA5007,4959 emission at z = 2.393. 

MACS J2135-0102 For completeness, we also include in 
our sample discussion of the velocity field of the lensed Ly- 
man Break Galaxy (LBG) at z = 3.07 in MACS J2135-0102 
(also kn own as the "Cosm ic Eye" due to its lensed mor- 
phology; IStark et al.l l2008h . Multi- wavelength studies have 
shown that the source is a L* Lyman-break galaxy with a 
stellar age of 80-300 Myr, star formation rate of SFR = 40- 
TOM^yr" 1 , and st e llar mass of M* = 6 ± 2 x 10 9 M 
ICoppin et al.ll2007i ; ISiana et alj|2009h . 

2.2 OSIRIS Observations and Data Reduction 

Detailed 2-D spectroscopic observations of the targets listed 
in Table 1 were made with the near-infr ared integral field 
spectrograph OSIRIS l|Larkin et al.ll2006h on the 10 meter 
Keck II telescope usin g the laser guide star a daptive op- 
tics (LGSAO) system 1 Wizinowich et all [2006) to correct 
for atmospheric distortion. A suitably bright star (_R < 17) 
within ~ 50" of the target was used for tip-tilt correction. 
We used the 100 milli-arcsecond pixel scale in all observa- 
tions which provides a field of view of at least 3.2" x 6.4". 
Narrow band H and K filters (Table 1) were selected to 
target the Ha or [Oiii]A5007 emission line at a spectral res- 
olution R ~ 3600, corresponding to FWHM ~ 6 A in the K 
band. MACS J2135-0102 and MACS J0451+0006 were ob- 
served on 2007 September 02 and 03 whilst the remaining 
targets were observed between 2008 November 27 - 30. 
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All observations were taken in 0.6 — 1.3" seeing. Before 
observing each arc, we took short exposures of the tip-tilt 
reference star to center the IFU pointing. The tip-tilt ex- 
posures were also used for flux calibration and to calculate 
the point spread function. Gaussian fits to the point spread 
functions of the tip-tilt stars yield 0.13-0.20" FWHM res- 
olution for the 100 milliarcsecond pixel scale. With smaller 
pixel scales, the resolution delivered by the LGSAO system 
was 0.11" FWHM. Observations of each target were done 
in a standard ABBA position sequence to achieve good sky 
subtraction. In the case of C10024+1709 we chopped 8" to 
sky, whilst in the remaining cases we chopped the galaxy 
within the IFU. Individual exposures were 600-900 seconds 
and each observing block was 2.4-3.6 ks which was typically 
repeated three to six times. The total integration time for 
each object is given in Table 1. 

Our data redu ction methods clos ely followed those de- 
scribed in detail bv lStark etafl (120081). We used the OSIRIS 
data reduction pipeline ( Larkin et al.l|20 06) to perform sky 
subtraction, spectral extraction, wavelength calibration, and 
form the data cube. To accurately combine the individual 
data cubes, we created images of the integrated emission line 
and used the peak intensity to centroid the object. We then 
spatially aligned and co-added the individual data cubes to 
create the final mosaic. Flux calibration was performed by 
equating the flux density of the tip-tilt stars measured from 
2MASS photometry with the observed OSIRIS spectra. We 
estimate that the uncertainty in flux calibration is typically 
10%. 



2.3 Gravitational Lens Modeling 

In order to investigate the source plane properties we must 
first correct for the distortion and magnification by the 
cluster lens. We summarize here the ingredients necessary 
to construct the relevant cluster mass models. We will 
follo w the methodology def i ned by earlier re l evant arti- 
cles (iKneib et all Il993l . Il996l ; ISmith et all 120051 ; I Julio et all 
l2007h within which further details can be found. 

Our bas i c ap p roach is to use the code LenstoolQ 
l|Kneib et al.l Il993l ; lJullo et all 120071 ) to constrain a pa- 
rameterized model of the dark matter distribution. For 
each cluster, the model comprises two components: 1 or 
2 cluster-scale dark matter halos, parametrized with a 
dual pseudo-isotherma l elliptical mass distribution (dPIE, 
lElfasdottir et all 120071 ). and ~50 galaxy-scale dPIE dark 
matter halos, centered on massive cluster members occu- 
pying the strong lensing region in order to account for the 
presence of substructure. These galaxy-scale halos are as- 
sumed to have mass properties that follow a scaling relation 
based on the luminosity of the und erlying galaxy, ass uming 
a constant mass-to-light ratio (e.g. ISmith et al]|2005l) . 

Details on the c onstruction of eac h gravitational lens 
mode l are given in [5ve et all (|2007l ) for MACS J2135- 
0102, [Limousin etafl (l2009h fo r C100 24+1709 and MACS 
J0744+3927 and iRichard et alj i|2009] ) for the remaining 3 



sources in Table 1. Strong lensing constraints originate from 
the identification of 2-5 multiply-imaged systems per clus- 
ter within the various ACS images. We use the astrometric 
positions and spectroscopic/photometric redshifts of these 
sourcefl as individual constraints to derive best fit parame- 
ters on the mass distribution. Lenstool uses a Markov Chain 
Monte Carlo (MCMC) sampler to derive a family of mass 
models suitably fitting the strong lensing constraints, and 
we use these to derive the uncertainty on each parameter of 
the mass distribution. 

For each of the sources presented in this paper, the best 
model was then used to derive the geometrical transforma- 
tion necessary for mapping the source plane coordinates into 
the image plane. This transformation enables us to recon- 
struct the HST morphology and Ha emission line images 
in the source plane assuming conservation of surface bright- 
ness. The spot magnification fi xy and its associated error are 
computed with Lenstool at different positions across the 
object, using the family of mass models from the MCMC 
sampler. We can verify this value by computing the total 
magnification from the ratio of the sizes (or equivalently, the 
total fluxes) between the image and its source plane recon- 
struction. As the magnification factor is not isotropic, the 
angular size of each image is more highly stretched along a 
specific orientation (Figure[TJ thus affecting our source plane 
resolution. The linear factors /ii and \ii of the magnification 
(with fj, = fii x iai) together with their associated errors are 
listed in Table 1. A detailed illustration of the uncertainties 
of the mass modeling method for MACS J2135-0102 is given 
in lStark et all (|2008l i. 

A key parameter in our analysis is the physical reso- 
lution we achieve in the source plane for each target. To 
measure this, we use observations of the tip-tilt reference 
stars. These serve this purpose well as they are point sources 
observed with conditions and an instrumental configuration 
identical to those of our distant targets. We "reconstruct" 
these stars in the source plane as if they were located at the 
arc position using the same transformation as for the lensed 
galaxies, and fit a bivariate Gaussian to the point spread 
function in order to determine the source plane resolution. 
The typical FWHM of each arc in the direction of highest 
magnification is listed in Table 1 and, with the exception 
of C10024+1709 which is not highly-magnified, varies from 
60-350 pc with a mean of 200 pc. 



3 ANALYSIS 

First, we reconstruct the reduced, flux-calibrated data cubes 
to the source plane using transformations from the grav- 
itational lens models described in §2.3. The resulting data 
cubes were binned such that each spatial pixel corresponded 
to 0.5-1 FWHM resolution elements along the direction of 
highest magnification (Table 1). We fit Gaussian profiles to 
the strongest emission line (Ha or [Oiii]) at each spatial 



1 http : //www. oamp. f r/cosmology/lenstool 



Spectroscopic redshifts are available for a majority of the sys- 
tems 
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Table 1. 
Target List 



Name 




Q J2000 


<5j2000 


z 


texp 


Emission lines 


Mi 


X 


/'2 




FWHM 


HST ACS 






h m s 


o 




(ks) 


and OSIRIS filter 










(PC) 


photometry 


CI 0024+1709 




00 26 34.43 


+1709 55.4 


1.680 


16.5 


Ha, [Nil]; Hn5 


0.8 


X 


1.7 


1.38 ±.15 


820 


B,V,r,i,z 


MACS J04514 


-0006 


04 5157.27 


+00 06 20.7 


2.014 


14.4 


Ha, [Nil]; Knl 


1.3 


X 


37.0 


49 + 11 


60 


V,I 


MACS J07124 


-5932 


071217.51 


+59 3216.3 


2.648 


16.2 


Ha, [Nil]; Kc5 


1.5 


X 


18.7 


28 + 8 


90 


v.i 


MACS J0744-t 


-3927 


07 44 47.82 


+39 27 25.7 


2.209 


14.4 


Ha, [Nil]; Kn2 


1.9 


X 


8.6 


16 + 3 


310 


v.i 


C10949+5153 




09 52 49.78 


+5152 43.7 


2.394 


19.2 


[Om]; Hn4 


1.2 


X 


6.0 


7.3 + 2.0 


350 


V 


MACS J2135-0102 


213512.73 


-010143.0 


3.074 


21.6 


H/3, [Om]; Knl 


3.5 


X 


8.0 


28 + 3 


120 


V.I 



Table 1. Notes: Gravitational lens modeling is discussed in Section 2.3. The lensing amplification is non-uniform and highly directional, 
thus we give the typical linear magnifications fii, 112 along the minor and major lensing axes as well as the overall flux magnification 
fi. The source plane resolution refers to the typical FWHM of a point source in the direction of highest magnification. The photometric 
bands B, V, r, i, I, z are ACS filters F435W, F555W, F625W, F775W, F814W, and F850W respectively. 




MACSJ0712 + 5932 z = 2.64 



3 2 1 0-1-2-3 




CI0949 + 5152 z = 2.38 



t 

t 
f 



2 0-2 




Figure 1. False color ACS images of the sample derived from the multi-band data referred to in Table 1. The xy scale is in arcseconds. 



pixel using a weighted % 2 minimization procedure and de- 
termine the two-dimensional intensity, velocity, and velocity 
dispersion maps. To compute the emission line fits we first 
subtracted the median value at each spatial pixel to remove 
any source continuum and residual sky background. A blank 



region of sky within each data cube was used to determine 
the sky variance spectrum V(X). The spectra to be fit are 
weighted by w(X) = V -1 appropriate for Gaussian noise, so 
that regions of higher noise (e.g. strong sky emission lines) 
do not cause spurious fits. We compute the % 2 statistic for 
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the best-fit Gaussian as well as a for a featureless spectrum 
(/(A) = 0) and require a minimum improvement over the fit 
with no line of A\ 2 = 16-25 for the various arcs (i.e. 4-5<r 
emission line detection). If this criterion was not met, we av- 
eraged the surrounding 3x3 spatial pixels to achieve higher 
signal-to-noise. No fit was made if the 3x3 averaging still 
failed to produce the minimum A\ 2 improvement. We cal- 
culated the formal la error bounds by perturbing the Gaus- 
sian fit parameters until the \ 2 increases by 1 from the best- 
fit value. In all following sections, we have deconvolved the 
line widths with the instrumental resolution (R ~ 3600) by 
subtracting the instrumental resolution in quadrature from 
the best-fit Gaussian a. The resulting source plane inten- 
sity, velocity, and dispersion fields for the entire sample are 
shown in Figure [5] and demonstrate detailed kinematic and 
morphological properties on scales down to 100-200 pc. 

From the source plane nebular emission line intensity 
and dynamics we estimate the size, dispersion and dynami- 
cal mass of each galaxy. The size is calculated as the max- 
imum diameter from pixels with successful emission line 
fits, roughly equivalent to a major axis diameter with limit- 
ing isophote ~ 10~ 16 ergs" 1 cm" 2 arcsec" 2 . The uncertainty 
is dominated by errors in the lensing magnification and is 
< 20%. We estimate the global average velocity dispersion 
of each galaxy as the flux-weighted mean of the fit pixels, 
(7 mM „ = ~}2 Vpixlpix with typical uncertainty of ~ 5% esti- 
mated from errors in the fit parameters and the spread of 
a in individual pixels. This measurement of the dispersion 
is not affected by resolved velocity gradients. The six galax- 
ies in our sample all have large (T me „„ = 50 — 100 kms" 1 , 
which is consistent with other resolved observations of non- 
lensed galaxies at z £ 2 (e.g. Forster Schreiber et al.l 120091 : 
lLaw et al.ll2009l ; iLehnert et alj|2009h 



3.1 Kinematic Modeling 

To test whether the kinematics of each galaxy are consistent 
with a rotating system, and to estimate the inclination of 
any disks, we construct and fit simple disk models to the 
observed velocity fields. The d isk model for MA CS J2135- 
0102 is discussed in detail by IStark et al.l (|2008l ). and we 
follow a similar method for the rest of the sample. We use 
an arctangent function to estimate the circular velocity as a 
function of radius, 



V(R) = Vo + - V c arctan ^ 

7T Rt 



(1) 



which ICourteaul (|l997h showed to be an adequate simple 
fit to galaxy rotation curves in the local universe. The disk 
models contain seven parameters: inclination i, position an- 
gle 9, coordinates (a, 5) of the disk center, scale radius R t , 
asymptotic velocity V c , and systemic velocity Vo- 

To test how well these simple models can describe the 
data, we constructed velocity fields covering a large range 
of parameter space for each galaxy. From the disk center, 
inclination, and position angle we compute the disk radius 
at each pixel in the source-plane velocity maps. We then cal- 
culate the circular velocity from Equation [T] and correct for 
the azimuthal angle and inclination to obtain the observed 



velocity field of the model disk. To simulate lensing distor- 
tion, the models were then convolved with an elliptical point 
spread function estimated as the best 2D Gaussian fit to the 
reconstructed tip-tilt reference star images. To estimate the 
goodness-of-fit we compute the \ 2 statistic using the la er- 
ror from emission line fits. The emission line fits routinely 
yield a reduced x 2 close to unity indicating that these errors 
are a good estimate of the true uncertainty. To estimate un- 
certainty in the parameters we perturb the model until the 
X 2 increases by one standard deviation from the best fit. 
Best-fit \ 2 values for each galaxy are given in Table 2, ex- 
cept for CI 0949+5152 for which no reasonable fit was found. 
Velocity contours for the best-fit disk models are shown in 
Figured 

For galaxies which can be reasonably well described 
with disk-like kinematic structure, we extract the rotation 
velocity, inclination and disk position angle (major axis) 
from the model. The latter two values are usually estimated 
from morphology, but such methods are not necessarily ac- 
curate in this case due to the lensing distortion and asym- 
metric light distribution of these objects. We therefore adopt 
the inclination and position angle from best-fit disk models. 
In Figure[3]we show the one-dimensional velocity and disper- 
sion profiles extracted from the data along a slit aligned with 
the disk position angle (or with the morphological major axis 
for the merger CI 0949+5153). This allows us to better quan- 
tify the kinematic structure and, particularly, to estimate an 
approximate rotational speed Vmax for each galaxy. The cir- 
cular velocity is a parameter in the disk models, however, 
in some objects the velocity field does not reach an asymp- 
totic value (e.g. MACS J0744+3927) . In such cases the best 
fit value of V c does not necessarily accurately represent the 
observed range of velocity. However, 4/5 galaxies in addi- 
tion to MACS J2135-0102 display a coherent and roughly 
monotonic velocity curve consistent with the disk model. 
From these we estimate V ma x sin i as half of the range in 
the one- dimensional model profile in the region detected in 
Ha emission, which provides a more reliable measure of the 
observed velocity range than the best-fit V c sini. Values of 
Vmax sin i as well as the best-fit inclination from disk models 
are given in Table 2. 

To estimate the dynamical mass for each galaxy, we 
adopt two approaches. First, we use the velocity disper- 
sion and maximum extent of the Ha emission and calcu- 
late M a y n (R) = C Ra mean /G with C = 5 appropriate 
for a sphere of uniform density and the radius is taken as 
half the maximum diameter. The geometrical factor C is 
likely between 3 and 10 for these highly turbulent galaxies 
|Erb e t al. 2006) . The dynamical mass is likely uncertain to 
within a factor 2x due to the assumed geometry and ex- 
clusion of rotationally-supported mass, whereas random un- 
certainty propagated from R and owan is < 25%. We also 
use the velocity gradient to estimate the dynamical mass 
via Mdyn sin 2 i = R (V max sin i) 2 /G but note that the max- 
imum velocity shear in the sample gives A^/2(T mm „ g 1.5 
for all galaxies, so the rotational mass is less than half of the 
dispersion-supported mass before correcting for inclination. 
In future sections, we therefore typically adopt the dynam- 
ical mass calculated from velocity dispersion, noting that 
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the total uncertainty is roughly a factor of two dominated 
by the assumed C = 5 with an additional systematic error 
due to the ignored rotational support. Adopted values of 
the diameter, dispersion, and dynamical mass are given in 
Table 2. 

Before we discuss what can be learned from the sample 
as a whole, we briefly review the results for each of the six 
galaxies in turn. 



3.2 Review of Individual Properties 

8.2.1 CI 0024+1709 

After accounting for lensing amplification, the extent of Ha 
emission is ~ 10 x 18 kpc from our lens model, making 
this the largest object in our sample. The integrated Ha 
emission line flux implies an intrinsic SFR= 27 ± 6 M©yr -1 
(uncorrected for dust extinction). 

The Ha emission line and rest-UV morphologies both 
show a clumpy ring surrounding a fainter central region 
(Figure [2]), similar to that seen in some other massive, ex- 
tended galaxies a t high redshift (e.g . the z = 2.26 galaxy 
Q2643-B X482 in ICenzel et all [20081 ) as well as local ring 
galaxies l|Romano et al. 2008T I. The narrow ring, ~10kpc 
diameter, and misaligned central feature suggest that the 
ring of concentrated star formation was triggered by an on- 
axis collision with another galaxy, as opposed to dynamical 
resonances within the galaxy. Thi s hypothesis is supported 
by the numerical simulations of iLvnds fc Toomref (|l976l ) 
and id entification of i nterac ting companions to ring galax- 
ies by iRomano et al.l (|2008T ) . Unfortunately, the Ha emis- 
sion in the central region coincides with a detector glitch 
in most of the IFU exposures and so the central disper- 
sion and star formation density are only poorly constrained, 
with a = 50 ± 21 kms -1 . In the ring, five UV- and Ha- 
bright knots are resolved with similar isophotal radius, ve- 
locity dispersion and luminosity of R ~ 0.5 kpc (for Lhc > 
1.5 x 10 -16 ergs -1 cm -2 arcsec -2 ), a ~ 45 — 65 kms" 1 , and 
Lhc ~ 2 x 10" 17 ergs" 1 cm' 2 , with uncertainties from the 
lens model and emission line fits of ~ 40% in radius, 15% 
in dispersion and 30% in luminosity. This implies a dynam- 
ical mass of the Ha clumps from the velocity dispersion of 
M dyn = 2±lxl0 9 M assuming uniform density, and a Ha 
star-formation rate of SFR = 2.0 ± 0.7 M©yr -1 . The veloc- 
ity field shows a smooth velocity gradient around the ring 
with peak-to-peak amplitude 2 V c sini = 120 ± 10 kms -1 . 
Using the axial ratio from the optical imaging, we estimate 
cosi = 0.6 ±0.1 which suggests a rotationally supported 
mass M dyn = RV^/G = 2.1 ± 0.8 x 10 10 M© within a ra- 
dius of 9 kpc. If the galaxy has a uniform spherical mass 
distribution, then the observed dispersion suggests a larger 
dynamical mass of 5 x 10 10 Mq. The resolved star form- 
ing regions in the ring thus account for about 50% of the 
integrated Ha flux and 15-25% of the mass interior to the 
ring. 



3.2.2 MACS 30451+0006 

The integrated Ha emission line flux suggests a star- 
formation rate of SFR=7 ± 2 M©yr -1 , with the two bright- 
est clumps (Figure [2} contributing 0.8 ± 0.2 and 3.6 ± 
0.9 M©yr -1 . Some parts of the galaxy were not observed 
with OSIRIS due to the limited field of view and 20 arcsec- 
ond extent of the arc. We expect these regions to contribute 
an additional 3 ± 1 MQyr" 1 based on their rest-frame UV 
flux. The large magnification of this galaxy (/i = 49 ± 11) 
leads to an unprecedented 60 pc typical resolution in the 
source plane, providing unique insight into the kinematic 
structure of the star forming regions. The Ha velocity map 
reveals a coherent velocity field across the galaxy and within 
the clumps with a peak to peak amplitude of 80 ±20 kms -1 
(Figure|3]). The bi-symmetric velocity field together with the 
good disk model fit (x 2 =3.6) suggest that the kinematics of 
this system are likely dominated by rotation rather than 
merging. 



3.2.3 MACS 30112+5932 

From our IFU observations, only the brightest UV clump is 
detected in Ha with an intrinsic star-formation rate of SFR 
= 5±1 M©yr -1 and we place limits on the remaining clumps 
of SFR < 3 M©yr -1 (3cr). The Ha flux in the brightest 
clump is consistent with the total galaxy flux derived from 
longslit spectroscopy, implying very little star formation in 
the rest of the system. The two prominent UV-bright regions 
have FWHM sizes <, 0.5 kpc and are separated by 1.5 kpc 
(Figure [2}. The Ha emission in the brighter region exhibits 
a velocity gradient of 38 ±8 kms -1 along the highly sheared 
direction. While the disk model provides an excellent fit to 
the velocity field with xi = 0.4, it is unclear whether the 
dynamical state of MACS J0712±5932 is dominated by ro- 
tation, merging, or dispersion since only the strongest Ha- 
emitting region is detected with OSIRIS. Assuming a uni- 
form density within a radius of 0.2 kpc and a = 80±5 km s -1 
we derive a dynamical mass of M dyn = 1.4 ± 0.2 x 10 9 Mq 
for the nebular emission region. The N[ll]A6583 emission line 
is detected at the 2a level in individual spatial pixels con- 
sistent with the flux ratio N[ll]A6583/Ha = 0.2 observed in 
NIRSPEC longslit spectra, with no measurable gradient in 
the flux ratio. 



3.2.4 MACS 30144+3921 

Ha emission is detected across 2.0 kpc in the source plane 
and is resolved into 3 bright clumps in both the rest-frame 
UV continuum and Ha. The two dimensional velocity map 
clearly shows a bi-symmetric velocity field and is well de- 
scribed by the disk model with Xv = 3.8 and peak to 
peak amplitude of 280 ± 25 kms" 1 , the largest observed in 
any high-redshift galaxy of this size (Figure [3]). The ve- 
locity dispersion varies from ~ 60-130 kms" 1 and peaks 
near the dynamical center of the galaxy. This suggests a 
turbulent primitive disk with a total star-formation rate of 
2.4 ± 0.5 MQyr -1 . The dynamical mass within a radius of 
lkpc is M dyn =R V c 2 /G = 5 ± 1 x 10 9 M© /sin 2 i inferred 
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Figure 2. Source plane properties of the lensed galaxy sample. From left to right in each row: ACS broad- band emission, line intensity, 
velocity and velocity dispersion derived from the OSIRIS nebular emission line data. Axes are in kpc with 1-2 pixels per OSIRIS 
resolution element. Each galaxy shows morphological and kinematic structure on scales down to ~ 100 — 200 parsecs. The nebular line 
flux distribution is similar to the rcst-UV morphology with multiple resolved clumps, shown as red crosses on the ACS images. Red 
contours on the line intensity maps denote clumps which are spatially extended and unconfused at the ~ 3cr flux isophote in narrow-band 
emission line images, as described in §4. Best-fit disk model contours are shown on the velocity maps for those galaxies whose dynamics 
are consistent with rotation (§3.1). Dashed boxes in the ACS images of MACS J0451 and MACS J0712 indicate the smaller regions for 
which OSIRIS data are displayed. 



from the velocity gradient, or Md yn ~ 11 ± 2 x 10 M© 
from the flux-weighted mean dispersion assuming uniform 
density. 

The large rotation velocity, low SFR, and central dis- 
persion peak conceivably make MACS J0744+3927 the most 
evolved disk-like galaxy of its size yet observed at high red- 
shift and show that small, low mass (~ 10 10 M ), turbulent 



field disk galaxies are already in place only 3 Gyr after the 
big bang. 



3.2.5 C10949+5153 

The reconstructed image of CI 0949+5153 at z = 2.393 
shows two diffuse regions separated by ~ 3.9 kpc, with the 
southwest component redshifted by 17 ±6 kms -1 from its 
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companion (Figure [3}. Both components show resolved ve- 
locity gradients of order 30 and 100 kms" 1 but in opposite 
directions, indicating that these two regions are likely merg- 
ing (this system can not be described by a disk model). 
Each component is further resolved into multiple clumps 
of diameter 0.3 — 0.7 kpc within an isophote of Lyoiu] = 
1.8 x 10~ 16 ergs -1 cm~ 2 arcsec~ 2 . Interestingly, while the 
clumpy [Oni] morphology appears to correlate with the rest- 
frame UV, the ratio of rest-UV to nebular emission surface 
brightness varies considerably between clumps. This could 
be caused by non-uniform extinction or different star forma- 
tion histories, since the rest-UV traces older star formation 
episodes than Ha flux (;£ 100 and ;£ 20 Myr respectively). 
Assuming a global H/3/[Om] flux ratio from longslit spectra 



and case B recombination, the H/3-derived star formation 
rate of the entire system is 20 ± 6 Mq yr _1 with 63% of the 
nebular emission in the larger northeast region. 



Dynamically, it appears that CI 0949+5153 is the only 
major merger in our sample, containing two spatially sepa- 
rated components with opposite velocity gradients and en- 
hanced dispersion between components. From the velocity 
dispersion we crudely estimate that the northeast and south- 
west components have mass 12 x and 4 x 10 9 Mq respec- 
tively. Both of these components are fragmented further 
into multiple large clumps with seemingly coherent veloc- 
ity fields. 
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3.2.6 MACSJ2135-0102 

For completeness, we also include in our sample discussion 
of the velocity field of LBG J2135-0102 (also k nown as the 
"Cosm ic Eye" due to its lensed morphology). IStark et al.l 
(2008) show that the global dynamics suggest a thick galaxy 
disk in an early state of assembly with rotation speed (cor- 
rected for inclination) of V c — 67 ± 7 kms -1 and V c /a r = 
1.2 ±0.1. 

The source plane morphology of this galaxy comprises 
two distinct star-forming regions embedded in a rotating 
disk. The northern component has a diameter of ~ 0.8 kpc, 
a star formation rate of 21±3 Mq yr _1 and a peak line width 
of 50 ± 5 kms -1 . Together, these suggest an intensely star 
forming region with dynamical mass Mdyn = 1.2 ± 0.5 x 
10 9 Mq . The southern component is brighter in the rest- 
frame UV, but has a lower H/3-inferred star formation rate 
of SFR = 11 ±2 M yr _1 . The diameter of 1.3 ±0.3 kpc and 
peak line width of 65 kms -1 suggest a mass M — 3 ± 1 x 
1O 9 M . 



3.3 Properties of the Ensemble 

Having discussed the sources individually, we now turn to 
general conclusions based upon their integrated properties 
as listed in Table 2. 



It is illustrative to demonstrate the dramatic improve- 
ment in our understanding of the internal dynamics of our 
sources that arises uniquely through the improved spatial 
sampling enabled by studying strongly-lensed galaxies. To 
do this we bin the source-plane data cubes of our targets 
to a coarser resolution typical of that achievable for un- 
lensed sources observed with an AO-corrected resolution of 
0.15". We then re-fit the emission lines and extract one- 
dimensional velocity profiles using the same methods as 
adopted for the original data. Since the lensing magnifica- 
tion results in \x times more flux spread over /j, times as 
many pixels, the binned signal-to-noise ratio is a a factor 
yfji higher than equivalent observations of unlensed ver- 
sions of our galaxies. These simulations therefore repre- 
sent non-lensed galaxies observed for much longer integra- 
tions (~ 100 hours!). Even so, the resulting velocity pro- 
files (also shown in Figure [3j) are considerably inferior to 
those of our lensed data. A credible rotation curve is only 
retrieved for C10024±1709, our least magnified galaxy, with 
important kinematic detail lost in all other objects (for ex- 
ample, MACS J0712±5932 is unresolved). Such poor spa- 
tial sampling is insufficient to distinguish between rotation 
and merging. More quantitatively, velocity gradients in all 
of our objects except C10024±1709 are significantly under- 
estimated (typically Vmax,lens = 0.6 ±0.2). 



3. 3. 1 Kinematics 

All six objects in our sample show well-resolved veloc- 
ity structure. A visual inspection of Figures [5] and [3] re- 
veals that four galaxies (CI 0024±1709, MACS J0451±0006, 
MACS J0744+3927, and MACS J2135-0102) have smooth 
velocity gradients suggesting ordered rotation whilst 
C10949±5153 appears to be a merger with both com- 
ponents showing well-resolved velocity gradients. The dy- 
namical state of MACS J0712±5932 is unclear since only 
the brightest Ha-emitting region is detected with OSIRIS. 
The nebular emission shows a smooth velocity gradient in 
MACS J0712±5932, but subtends only four resolution ele- 
ments in the most highly magnified of the three images and 
so it is unclear wither the kinematics of MACS J0712±5932 
are dominated by rotation, dispersion, or merging. 

For the galaxies studied here, the best fit disk models 
have xt values which range from 0.4 to 5.9 with a mean of 
3.5, indicating agreement typically within 2a for individual 
pixels. For reference, disk models for 18 galaxies in the SINS 
survey which show the most promine nt rotation yi e ld bes t-fit 
xl values of 0.2-20 to the models of ICresci et al.l (|2009h . so 
our simple model provides an equally good fit to the lensed 
galaxies. All galaxies show small-scale deviations from the 
model as indicated by the typical Xu > 1; these proper mo- 
tions could be caused by the effects of gravitational insta- 
bility, or simply be due to the unrelaxed dynamical state 
indicated by high velocity dispersions a > 50 kms -1 . We 
therefore conclude that the model provides an adequate fit 
to the data and that the velocity fields are consistent with 
the kinematics of a turbulent rotating disk, except in the 
case of C10949±5152. 



3.3.2 Physical Characteristics: Size, Luminosity and Mass 

Next, we examine the integrated physical properties of the 
lensed galaxies. First we briefly compare the luminosity, size 
and mass to the general Lyman break population at z ~ 2- 
3. In Figure U we compare the distribution of apparent R- 
band magnitudes for our sample to that o f the Lyman-break 
population at z ~ 3 l|Shaplev et al.|[200ll ). The comparison 
demonstrates that five of our six galaxies are fainter than 
L* , ranging from 0.1 — 1.5 L* , with median 0.5 L* well be- 
low that of other spatially resolved IFU studies. The in- 
trinsic Ha flux of the arcs is also lower than in other sur- 
veys, with a mean and median inferred S FR of 17Mpyr -1 
(Ta ble 2) compared to 26-33 Mpy r -1 in lLaw et all (|2009f) 
and iForster Schreiber et al.l |2009). Comparing the typical 
spatial extent of the nebular emission to other surveys, 
the median radius (extent of detected nebular emission) 
and dynamical mass of the lensed sample are 2.1 kpc and 
1.3 x 10 10 Mq, some what more extend ed than the compact 
galaxies studied by lLaw et all (|2009h which have median 
radius of 1.3 kpc and M dyn = 0.7 x 10 10 M Q . Turning to 
the dynamics, the flux-weighted mean velocity dispersions 
o'mean are perfectly consistent with those studied by Law 
et al. (2009): both samples span a range of 50-100 kms -1 
with mean and median 70-75 kms -1 , indicating that the 
two sets of galaxies probe the same range of dynamical mass 
density. We note that the median Ha FWHM (from longslit 
spectra; lErb" et al.l [2006) and stellar mass of ~ 6 kpc and 
2.9 x 10 10 M are a better representation of the lLaw et al.l 
(2009) sample as they do not depend on sensitivity of the 
data. The SINS survey probes somewhat more extended ob- 
jects, with median Ha FWHM and stellar mass of ~8 kpc 
and M* ~ 3 x 10 10 Mq. The mass density and extent of neb- 
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Figure 3. Left to right: velocity field, one dimensional extracted velocity profile, and equivalent velocity dispersion profile. Velocity 
and dispersion curves are extracted from left to right along the slit shown on the left figures. The rotation curve one would obtain in 
the absence of lensing with 0.f5" resolution typical of AO-corrected observations is shown by the blue diamonds. Best-fit profiles of the 
velocity field are shown in red. No model is shown for Cl 0949+5153 since the velocity field is inconsistent with a rotating disk. 
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Table 2. 
Intrinsic Galaxy Properties 



Galaxy 


Diameter 


L/L* 


SFR 


Mdyn 


(7 


Vmax sini 


i 


X.7j 




(kpc) 




(Moyr- 1 ) 


(10 9 M S ) 


(kms- 1 ) 


( km s — 1 ) 


n 




CI 0024+1709 


20 ± 2 


0.7 


27±6 


55±10 


69±5 


76±12 


50+ 10 
ou -30 


5.9 


MACS J0451+0006 


5 ± 1 


0.4 


7±2 


19±4 


80±5 


38±10 


9^+20 
zo -10 


3.6 


MACS J0712+5932 


1.5 ±0.3 


0.2 


5±1 


5.7±1.2 


81±2 


15±4 




0.4 


MACS J0744+3927 


2.0 ±0.3 


0.1 


2.4±0.5 


11±2 


99±4 


129±12 




3.8 


CI 0949+5152 


7 ± 1 


0.5 


20±6 




66±3 








NE 


4.0 ±0.6 




13±5 


12±2 


71±2 








SW 


2.0 ±0.3 




7±3 


3.8±0.6 


57±2 








MACS J2135-0102 


3.5 ±0.3 


1.5 


40±5 


5.9±1.0 


54±4 


47±7 


55 ± 10 


3.7 



Table 2. Source-plane properties of the lensed galaxies. Values for MACS J2135-0102 are reported in Stark et al. (2008). Diameter is the 
extent of nebular emission except for MACS J0712+5932 where it represents the extent of detected V-band emission, with uncertainty 
accounting for spatial resolution and errors in the lens model. The luminosity is compared to L* for the LBG population at z = 3 (see 
text for methodology). Error in luminosity is 15—30% dominated by the uncertain magnification. Velocity dispersion a is taken as the 
flux-weighted mean. Vmax sin i is computed as half the maximum peak-to-peak velocity from fits to the velocity field. Inclination i is the 
value relative to the plane of the sky for the best-fit disk model, with uncertainty given as the range of parameter space within 1<t of the 
best fit. CI 0949+5152 is separated into the northeast and southwest merging components; no circular velocity is given since rotating disk 
models do not provide a good fit. Dynamical mass is calculated as C Rcr 2 /G using the diameter and dispersion values given in the table 
with C = 5 and error from R and cr, although the true uncertainty in M ayn is roughly a factor of two due to unknown mass distribution 
and rotation support. 




Figure 4. Apparent R-band luminosity function of z ~ 3 
LBGs from Shapley et al. (2001) compared to the distribution 
of suitably-corrected apparent R magnitudes for our lensed sam- 
ple as well as the SINS (Forster Schrcibcr et al. 2009) and Law 
et al. (2009) samples. We compute the apparent R magnitudes of 
the IFU samples at z = 3 for an effective wavelength A = 600—800 
nm depending on the available photometry. 



ular emission of the len sed galaxie s are th erefore comparable 
to the more luminous lLaw et all l|2009l ) sample and some- 
what lower than in the SINS survey. The lensed galaxies are 
also below the average size and mass for L* systems at sim- 
ilar redshifts, with the notable exception of CI 0024+1709. 

Our high resolution data enable us to determine re- 
liable dynamical properties and compare with other IFU 
observations. In particular we can examine the preva- 
lence of dispersion a nd rotation as a function of size 
and dynamical mass. lLaw et all l|2009h report evidence 
for dispersion-dominated kinematics in compact low-mass 



galaxies in contrast to the ~equal mix of rotation, disper- 
sion, Mid_£omrjle2£/merg^r_jdnernatics found in the SINS sur- 
vey (|Forster Schreiber et al.ll2009h , although the dispersion- 
dominated fraction is higher for more compact SINS galax- 
ies. However, this claim is hampered by the small number 
of resolutio n elements (^2-4) s ubtended by each source (al- 
though see lEpinat et al.ll2009T ). We can address this issue 
with the superior resolution offered by gravitational lens- 
ing. The relevant relations between V/o and size or dynam- 
ical mass is shown in Figure [S] for the five lensed galaxies 
which are reasonably w ell fit by d i sk mo dels compared to 
the compact LBGs from lLaw et alj d2009f) and the rotation - 
dominated SINS galaxies described bv lcresci et ail l|2009[h 
showing that the ratio of velocity shear to dispersion is 
higher for larger, more massive galaxies. Furthermore, the 
general agreement between the lensed sample and other 
comparison samples is also shown in Figure [5] and demon- 
strates that the dynamics of sub-L* star-forming galaxies do 
not differ substantially from more luminous objects. 

Noting that many galaxies do not reach an asymp- 
totic velocity within the region probed by our IFU obser- 
vations, it is likely that larger circular velocities exist in 
faint outer regions of these galaxies. Indeed, the spatial ex- 
tent and mass inferred from longslit spectra and photom- 
etry are higher than from relatively shallow IFU data. We 
therefore also compare the observed velocity gradient and, 
where appropriate, inclination-corrected rotational velocity 
as a function of size and dynamical mass in Figure [5] The 
observed velocity gradients of 3/5 lensed galaxies and the 
entire comparison sample are within 0.3 dex of the median 
value Vshear/R = 25 km s~ 1 kpc - 1 . The consistent velocity 
gradients and corresponding V s hear/c obtained with a typ- 
ical a — 75 kms -1 (Figure O suggest that the extended 
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structure in some of the compact objects may in fact host 
velocity fields comparable to the larger disk galaxies. 

While MACS J0451+0006 and MACS J2135-0102 have 
velocity gradients consistent with the other samples con- 
sidered here, the clear outlier MACS J0744+3927 demon- 
strates that at least some small galaxies at this redshift 
have significant dynamical support from ordered rotation. 
This galaxy has an inclination-corrected V c /a = 1.8to'4 at 
radius R = 1 kpc and a Ha-derived star-formation rate sur- 
face density Y*sfr = 0.8 ± 0.3, which is a factor three lower 
than any galaxy in the dispersion-dominated sample, possi- 
bly indicating a later stage of evolution. 

Motivated by the relatively large velocity shear and 
small T, S fr of MACS J0744+3927, we explore the possibil- 
ity that the low V/a observed in many compact galaxies 
is an effect of selecting galaxies in early stages of evolution 
with large dispersions and star formation rates. This is ex- 
pected in models of star formation feedback, which predict 
that mechanical energy introduced by supernovae leads to 
higher velocity dispersion an d hence lower V/ a as T,sfr in- 
creases. In one such model, iDib et all (I2006T ) simulate the 
velocity dispersion within the inter-stellar medium of galax- 
ies induced by supernovae energy injection (assuming a 25% 
feedback efficiency) and predict a correlation between star- 
formation surface density and velocity dispersion. 

To test whether supernova feedback might be responsi- 
ble for the high dispersion observed in all z £ 2 star-forming 
galaxies, we show the velocity dispersion a as a function of 
Y^sfr in Figure[6] For the range of T>sfr with modeled su- 
pernova feedback, the observed velocity dispersions a = 30- 
80kms _1 are much higher than the predicted a = 10- 
15kms _1 . We note that the density used in these simula- 
tions is ~ 30 x lower than inferred for the z ~ 2 galaxies, so 
the simulated dispersion is likely overestimated by a factor 
of ~ 5. While the data show a slight trend of a increasing 
with Esfh, the observed relation is inconsistent with sim- 
ulations. We therefore conclude that supernova feedback is 
insufficient to explain the observed velocity dispersions. Fig- 
ure[6]shows the relation between V / o-Yjsfr, clearly demon- 
strating that V/a decreases with Y>sfr- However, this trend 
is likely ultimately due to the velocity-size relation: larger 
galaxies tend to have larger rotation velocities (Figure[5} and 
smaller Y>sfr (Figure [HJ. This is explained by different sen- 
sitivities of the data, as deeper spectra (lower T^sfr) reveal 
more extended structures at larger radius with larger rota- 
tion speed. The velocity-size correlation contributes much 
more to the observed V/a-T,sFR relation than any correla- 
tion between 'Esfr and the velocity dispersion. These data 
thus do not support the hypothesis that V/a is strongly af- 
fected by the density of star formation: a increases by less 
than a factor of 2 over two orders of magnitude in Y,sfr- 



Evidence is accumulating that the mode of star forma- 
tion may be very different in early s ystems compared to that 
seen locally l|Bournaud et alj [2008 : El megreen fc Elmegreenl 
2005). Rather than forming stars within giant molecular 
clouds which condense out of a stable galaxy, star forma- 
tion may be triggered by fragmentation of a dynamically 
unstable system. Briefly, in a rotating disk of gas and stars, 
perturbations smaller than a critical wavelength L max are 
stabilized against the inward pull of gravity by velocity dis- 
persion while those larger than some L m i n are stabilized by 
centrifugal force. If the dispersion and rotation velocity are 
too low, L m in > L max and perturbations of intermediate 
wavelength grow exponentially. This interplay is summa- 
rized by the Toomre parameter Q = L ma x/ L m i n which is 
calculated from th e velocity disp ersion, rotation curve, and 
mass distribution |Toomre Il964h . Galaxies with Q < 1 are 
therefore unstable on scales between L max and Lmin and 
will fragment into giant dense clumps. This could trigger 
star formation in clouds of much higher mass and radius 
than GMCs in local spiral galaxies with Q > 1, and can 
explain the clump-cluster and chain morphologies observed 
in many high-redshift galaxies. Dynamical friction, viscosity 
and tidal interactions may cause these clumps to migrate 
toward the center of the galaxy potential, forming a bulge 
which stabilizes the system against further fragmentation. 

From the galaxies whose velocity fields can be reason- 
ably well described by rotating systems, we calculate the 
Toomre parameter via: 



Q = 



ttGE 



(2) 



which describes the stability of a rotating disk of gas. If 
Q < 1 the system is unstable to local gravitational collapse 
and will fragment into overdense clumps. The value of k 
is somewhat uncertain as it depends on the unknown mass 
distribution; our observations are consistent with a range 
\/2-%— 2-^jr corresponding to constant V c and V c oc R respec- 
tively. Adopting k = V3Vmax/R appropriate for a uniform 
disk and using dynamical mass to estimate the surface mass 
density E, we find an inclination-corrected Q ^ 0.6 for all 
galaxies in our sample. We estimate that the uncertainty 
in Q is dominated by a factor of ~ 2 error in the dynam- 
ical mass. The assumed k introduces a negligible 15% un- 
certainty, with an additional random error of ~ 30% from 
the input parameters. Disk thickness and stellar abundance 
also affect the value of Q. Combined, these effects result in 
roughly a factor of 2 uncertainty. Even so, these galaxies all 
appear to be dynamically unstable since Q < 1. Hence we 
expect them to fragment into massive clumps on scales of 
order the Jeans length for dispersion support. In a uniform 
disk, the largest scale for which velocity dispersion stabilizes 
against gravitational collapse is 



3.4 Star-Formation Scales within Disks 

The small V ma x/a < 1.8 and clumpy morphology of all 
galaxies suggest that the rotating disks are highly turbu- 
lent and may be dynamically unstable. We therefore explore 
the scale lengths for gravitational collapse within the high 
redshift disk galaxies. 



Lj = 



8GE 



(•3) 



which is readily estimated from the dispersion and dynami- 
cal mass density. As with Q, the uncertainty in L j is a factor 
of ~ 2 dominated by the dynamical mass with additional un- 
certainty from the unknown mass distribution, disk height, 
stellar content, and directional dependence of a. The result- 
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Figure 5. Kinematics as a function of dynamical mass and ra- 
dius. Values for the lensed sample are given in Table 2 and de- 
scribed in the text. Data for the comparison samples are taken 
from Tables 3 and 4 of Law et al. (2009) and Table 2 of Cresci 
et al. (2009) and are discussed in £13.3,21 The median V sini/R = 
25kms~ 1 kpc -1 is plotted in the middle right panel, and the cor- 
responding V/cr for typical a = 75kms~ 1 is shown in the upper 
right. 



ing instability scale is 1-3 kpc for CI 0024+1709 and 0.1-1.5 
kpc for all other objects, consistent with the observed clump 
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Figure 6. Velocity dispersion, V s h ear / a , and radius as a function 
of star formation rate density estimated from Ho flux. Symbols 
are the same as in Figure [5] The solid line in the top panel is a fit 
to simulations; see text for details. The observed a is inconsistent 
with the simulations. Dispersion is at best weakly correlated with 
^SFR- The strong correlation of radius and cr/V s hear with T,$fr 
is likely due to different sensitivity of the observations. 



4 STAR FORMING REGIONS AT HIGH 
REDSHIFT 

For each galaxy in our sample, the gain in resolution pro- 
vided by gravitational lensing reveals multiple resolved giant 
Hn regions which, as a group, dominate the integrated star 
formation rate. This gain enables us to study, in more de- 
tail than hitherto, the properties of individual star forming 
complexes in typical high redshift galaxies. In the following 
analysis, we concentrate on those Hn regions with extended 
bright emission which are sufficiently isolated. Ideally we 
would construct a complete well-defined sample, however, 
in the exploratory work presented here we wish to limit the 
sample to objects for which the size and flux can be mea- 
sured robustly. We therefore specifically exclude Hn regions 
with overlapping emission in the Ha emission line maps 
(such as those in MACS J0744+3927). 

In brief, Hn regions were selected with a flux isophote 
in the image plane where the instrumental response is well 
known. Thus we can check whether each Hn region is suffi- 
ciently resolved by comparing with the point-spread function 
produced by a bright tip-tilt reference star. To measure the 



area, each Hn region was deconvolved with the point spread 
function using a method which accounts for the distorted 
shapes produced by gravitational lensing. The intrinsic flux 
and area were then computed by reconstructing the Hn re- 
gion in the source plane. 

Hn regions were selected above a flux isophote corre- 
sponding to ~ 3<r, where a is the noise level measured in 
narrow-band images, constructed by collapsing continuum- 
subtracted OSIRIS data cubes around the strongest emis- 
sion feature. Only isolated Hn regions (in which there is a 
single local maximum within the isophote) were selected. 
Regions not meeting this criterion are excluded from this 
analysis since attempting to disentangle multiple confused 
sources would introduce a large uncertainty. The isophote 
used varies from 1.4-4.6 x 10 -16 ergs -1 cm -2 arcsec -2 for 
different arcs. The isophote for each selected Hn region is 
shown in Figure [2] Each Hn region subtends 0.5-0.9 arcsec- 
onds in the image plane, with the notable exception of the 
giant clump in MACS J0451+0006 which spans 3.2 arcsec- 
onds. The selected Hn regions are all well-resolved since the 
point spread function has FWHM < 0.2" in all observations. 

For each Hn region we calculated the total nebular line 
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flux and intrinsic angular size in the image plane. The lu- 
minosity of each Hn region was taken as the total flux 
within the selection isophote, with the uncertainty deter- 
mined from the noise level in the narrow-band image. The 
Hn region size is more complicated to derive since the ef- 
fect of the spatial resolution needs to be taken into con- 
sideration. We proceed by making two basic assumptions. 
First, we assume the observed flux distribution is a con- 
volution of the actual distribution with the instrumental 
point spread function. We can therefore compute the ap- 
proximate intrinsic spatial variance by subtracting the Gaus- 
sian PSF variance from the measured variance within the 
isophote, Vint = Visa — &PSF- This is equivalent to subtract- 
ing a =FWHM/2.35 in quadrature from Gaussian profiles, 
but the procedure is applicable to non-Gaussian asymmetric 
distributions. In all cases Vi SO > &psf confirming that the 
Hn regions are well resolved: Vi so /(TpsF = 2.3 for one Hn 
region and > 3.2 for all others. Second, we assume that the 
isophotal area scales with spatial variance such that intrin- 
sic surface area A int — Aiso^^, where A iso is defined as 
the pixel size multiplied by the number of pixels within the 
isophote. Uncertainty in the isophotal flux ranges from 12- 
27% propagated from flux calibration and noise level in the 
narrow-band image. The uncertainty in A iso is estimated by 
assuming an error of PSF a for the isophote diameter. The 
range of resulting angular diameter (2y/ Ai so /w) and relative 
error is 0.12 ± 0.06 to 0.61 ± 0.08 arcseconds. 

To determine the source plane properties, each Hn re- 
gion is finally reconstructed using the lens model transfor- 
mation. The total magnification fi is calculated from the 
ratio of image- to source-plane flux. The source plane di- 
ameter is then defined as d — 2\ Aiso . This method is 

V A" 1 " 

robust against the spatially varying magnification of the 
lensed arcs. Uncertainty in the lens model has little effect 
on the diameter but significantly contributes to the error 
in flux. From the Ha flux (converted from [Om]A5007 for 
CI 0949+5 1 53) w e estimate the star formation rate using the 
iKennicuttl (1 19981 ) prescription. The values range from 0.8- 
3.6 M0yr _1 star formation rate and 0.3-1.0 kpc diameter. 

From the source plane diameters and kinematic data, 
we estimate the dynamical mass of the Hn regions. We use 
the flux-weighted mean velocity dispersion within the clump 
Omean an d assume a uniform spherical mass distribution 
(C = 5) as discussed in §3.1. The Hn regions have large flux- 
weighted mean dispersions, with a mean ranging from 45-80 
kms -1 . The dynamical masses range from 0.7-3 xlO 9 M©, 
far greater than typical star-forming clusters in the local uni- 
verse. Statistical uncertainty in the dynamical mass is 20- 
60%, so the additional factor of ~ 2 systematic uncertainty 
in C contributes significantly to the dynamical mass esti- 
mate. The mass-radius relation of the high-redshift clumps 
is shown in Figure[7]along with a sample of local star clusters 
and giant Hn regions for comparison. This shows that the 
Hn regions studied here are comparable in size and mass to 
the largest local star-forming complexes and consistent with 
the mass-radius relation observed locally. 

In Figure [8] we compare the size-luminosity relation of 
these star for ming regions with equivalent da ta observed 
locally (from iGonzalez Delgado fc Perezl 1 19971 for normal 



spirals, and iBastian et al.l 120061 for intense starbursts such 
as the Antennae) an d with a lensed arc at z = 4.92 
l|Swinbank et al.ll2009l ). The sizes observed in the high red- 
shift galaxies are comparable to the largest local star form- 
ing (Hn) regions but with ~100x higher star formation 
rates than in local spiral galaxies. However, the implied 
star formation rate densities are roughly consistent with 
that observed in the most vigorous local starbursts. This 
result is not significantly affected by the resolution of dif- 
ferent observations, although varying the selection isophote 
alters the radius-luminosity relation. We note that the flux 
isophotes used to define Hn regions in local spirals are much 
lower, so to quantify this effect we extract individual Hn re- 
gions from narrow-band Ha images of local galaxies from 
the SINGS an d 11HUGS surveys l|Kennicutt et al.l 120031 ; 
iLee et al.ll2007l ) using various isophotes. The effect on the 
size-luminosity relation as the isophote is increased is shown 
in Figure [S] as the vector A, which demonstrates that the ex- 
tracted regions become smaller with higher average surface 
brightness. The vector B shows the effect of degrading the 
resolution while the isophote is kept fixed, which causes dis- 
tinct Hn regions to merge together giving the appearance 
of a larger region with similar surface brightness. Since the 
surface brightness of the high-redshift Hn regions is ~ 100 x 
higher than in local spirals, we conclude that the offset can- 
not be explained by different resolution and sensitivity of the 
various data sets. Metallicity also fails to explain the offset, 
as demonstrated by MACS J2135-0102 which has an R23 in- 
dex and Civ P Cygni profile suggestin g 0.4-1 Z© metallicity 
l|Stark et al.ll2008l ; lQuider et al.ll2009h . This implies that the 
luminosities in the high redshift Hn regions are truly much 
larger than in local spirals. Furthermore the observed diam- 
eters are consistent with the Jeans length for support by 
velocity dispersion, suggesting that they collapsed as a re- 
sult of disk instabilities. We expect this to occur since the 
Toomre Q parameter is less than unity for all galaxies in our 
sample. It therefore seems likely that intense star-formation 
in high redshift galaxies is driven by the fragmentation of 
gravitationally unstable systems. 



5 CONCLUSIONS 

We have utilized strong gravitational lensing together with 
laser-assisted guide star adaptive optics at the Keck obser- 
vatory to study the internal characteristics of six z = 1.7 — 3 
star-forming galaxies with a spatial resolution down to ~ 
100 — 200 parsecs. Extending beyond the diffraction lim- 
ited capability of current 8-10 meter class telescopes, our 
study provides an interesting preview of the science that 
will routinely be possible with future 30 meter class opti- 
cal/infrared telescopes with adaptive optics systems (Ellis 
2009). We demonstrate that enhanced resolution allows us 
to resolve morphological and kinematic structure which has 
not been (and cannot be) discerned in AO-based surveys of 
similar non-lensed sources. In particular we resolve multi- 
ple giant star forming regions of size consistent with that 
expected from Toomre instability in all of our targets (in- 
cluding separate merging components). Furthermore we de- 
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Figure 7. Size versus dynamical mass of individual star-forming 
regions within our high-redshift sample, compared to local star 
clusters and Hn regions. The high-redshift Hn regions have size 
and mass consistent with the largest star-forming complexes in 
the local universe. Local data are taken from Harris (1996) and 
Pryor & Meylan (1993) for Galactic globular clusters; Hasegan 
et al. (2005), Hilker et al. (2007) and Evstigneeva et al. (2007) 
for ultra-compact dwarf galaxies; McCrady & Graham (2007) for 
M82 clusters; Fuentes-Masip et al. (2000) for giant Hn regions; 
and Bastian et al. (2006) for star-forming complexes in the An- 
tennae. The dashed line represents a model for regions which are 
optically thick to far-infrared radiation and have undergone adi- 
abatic expansion (Murray 2008). Dynamical mass for the high- 

T3 2 

redshift data is calculated as M,j.y n = C —q- with assumed C = 5 
appropriate for a uniform-density sphere. The mass error bars 
account for statistical uncertainty in R and <r. 



termine the dynamical state of galaxies which would be too 
poorly sampled to distinguish between rotation, merging, or 
dispersion-dominated kinematics without gravitational lens- 
ing. 

The lensed sample allows a unique study of relatively 
faint and small galaxies at high resolution. The median lu- 
minosity of the sample is half of the characteristic L* for 
z — 3 LBGs and well below that of other high-redshift IFU 
studies. We find that the kinematics of the sub-L* popu- 
lation are in general agreement with larger and more lu- 
minous galaxies at similar redshifts: four of the six sys- 
tems are clearly rotating with Vsini/a — 0.5-1.3, one 
has a velocity gradient consistent with rotation within the 
small Ha-bright region of the arc detected with OSIRIS, 
and one is a major merger. All have high velocity disper- 
sion a = 50-100 km s , consistent with all other resolved 
observations of z > 2 star-forming galaxies and compa- 
rable to the value a ~ 50km s~ for stars which forme d 
in the Galactic disk at z — 2 (|Sparke fc Gallagher! [2000). 
Galaxies with a larger radius of detected nebular emis- 
sion tend to have higher velocity shear and V/a and lower 
Ssfh; these trends are consistent with other observations. 
The case of MACS J0744+3927 with inclination-corrected 
V = 180^40° km s -1 clearly demonstrates that some small 
(R — TO ± 0.2 kpc), fast-rotating field disk galaxies are al- 
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Figure 8. Ha-inferred star formation rate as a function of Ho 
region size. Red circles represent individual Hn regions from our 
lensed sample. CI 0949+5152 was observed in [Oin] with the Ha 
flux taken to be 1.2 times the [Oni]A5007 flux as determined from 
long slit NIRSPEC observations. The value for MACS J2135-0102 
was determined from the H/3 flux assuming case B recombination. 
The star formation rates are all derived from Ha luminosities 
using the Kennicutt (1998) prescription. Red squares are from 
the lensed z = 4.92 galaxy presented by Swinbank et al. (2009), 
black points are Hn regions in local spiral galaxies from Gonzalcz- 
Delgado and Perez (1996), and blue open circles are star cluster 
complexes in the interacting galaxies NGC 4038/4039 from Bas- 
tian et al. (2006). The vectors A and B are explained in the text. 
The dashed line shows a fit to the data from local spiral galax- 
ies, and the dotted line is the same but increased by a factor of 
100 X in SFR density. The high-redshift star-forming regions have 
significantly higher density of massive stars than Hn regions in 
nearby spirals, but are comparable to the brightest star cluster 
complexes found in the local universe. 

ready in place by z = 2. Well-sampled velocity fields demon- 
strate that the ubiquitous giant clumps are often embed- 
ded in common rotating systems and are not independent 
merging components. In contrast, non- lensed clumpy sys- 
tems with velocity shear have been interpreted as mergers 
or remained ambiguous. These results demonstrate the need 
of high-resolution observations to distinguish rotating sys- 
tems and mergers, as well as to probe the < 1 kpc scale of 
star formation at high redshift. 

The observations presented here represent the best cur- 
rent probe of the scale of star formation at z > 2 and, as 
such, can be compared to theoretical predictions. Most of 
the star formation in our sample occurs within giant Hn 
regions of diameter 0.3-1 kpc comparable to the largest lo- 
cal Hn regions, with star formation rate density comparable 
to the most vigorous local starbursts and ~ 100 x higher 
than in typical spiral galaxies. This offset in T^sfr cannot 
be explained by the different resolution or sensitivity of low- 
and high-redshift observations. Only one of our six objects 
is a clear major merger with both merging components bro- 
ken into high-SFR clumps, suggesting that star formation 
episodes in sub-L* high redshift galaxies are triggered by 
Toomre instability independent of major merger events. We 
note that the merging system demonstrates no enhancement 
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in star formation (as traced by nebular emission) compared 
to the non-merging galaxies. The increased high-mass star 
formation relative to local spirals likely results from some 
combination of higher gas density, increased star formation 
efficiency, shorter star formation timescales, and possibly 
other effects. The clumpy star formation, high T^sfr, kine- 
matics consistent with rotation, and small Q < 1 are in 
qualitative a greem ent with the galaxy evolution model of 
iDekel et al.l <|2009h in which galaxies accrete their baryonic 
mass from cold streams. In this model the clumps form due 
to Toomre instability in turbulent disks and migrate into the 
galaxy center on timescales of a few hundred Myr, forming 
a bulge which stabilizes the system and increases the star 
formation timescale after ~ 2 Gyr. Numerical simulations 
suggest that cold-stream accretion is a dominant mechanism 
in galaxy assembly, especially at early times a nd in galax- 
ies w i th halo masses less than a few 10 11 M r l|Keres et al.l 
120091 : lOcvirk et alJliooi : iBrooks et all 120091). In part icular 
the high-resolution simulations of Brooks et al.l (|2009h show 



that accretion of cold gas dominates the buildup of stellar 
disks and bulges in sub-L* galaxies and allows disk forma- 
tion at earlier times than if the accreted gas is shock-heated. 
The rotating kinematics and distribution of star formation in 
the lensed sample therefore supports, at least qualitatively, 
a cold-stream accretion scenario for galaxy formation. 

Analysis of our sample has shown the benefit of using 
lensed sources to probe the inner structure of early star- 
forming galaxies. The advent of many panoramic multi-color 
imaging surveys will hopefully reveal many further examples 
of such systems, ensuring further progress with laser-assisted 
guide star adaptive optics. Ultimately, the next generation 
of 30- meter optical/near- infrared telescopes will extend the 
high-resolution studies presented herein to the entire popu- 
lation of high-redshift star forming galaxies. 
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